LPTHE-Orsay-96-78 



hep-ph/9610207 



October, 1996 



Power corrections in Drell-Yan production 
beyond the leading order 

G. P. Korchemsky f] 

Laboratoire de Physique Theorique et Hautes Energies f] 
Universite de Paris XI, Centre d'Orsay, bat. 211 
91405 Orsay Cedex, France 



Laboratory of Theoretical Physics 
Joint Institute for Nuclear Research 
141980 Dubna, Russia 



Abstract 

We study the power corrections to the cross-section of Drell-Yan production within the Wilson 
line approach and apply the methods of QCD asymptotic dynamics to identify the leading renormalon 
contribution. 



To appear in the Proceedings of the 
28th International Conference on High Energy Physics, 
Warsaw, Poland, 25-31 July 1996 



* E-mail: korchems@qcd.th.u-psud.fr 

t Laboratoire associe au Centre National de la Recherche Scientifique (URA D063) 



POWER CORRECTIONS IN DRELL-YAN PRODUCTION 
BEYOND THE LEADING ORDER 



G.P. KORCHEMSKY 

LPTHE, Universite de Paris XI, 91405 Orsay, France and LTP, JINR, 14-1980 Dubna, Russia 

We study the power corrections to the cross-section of Drell-Yan production within the Wilson line approach and 
apply the methods of QCD asymptotic dynamics to identify the leading renormalon contribution. 



1 Introduction 

To control the occuracy of perturbative QCD pre- 
dictions for physical observables it becomes im- 
portant to have a regular way of estimating cor- 
rections suppressed by powers of large momentum 
scale, Ap/Q p , and if not calculate explicitly the 
nonperturbative scale A p then find the leading ex- 
ponent p at least. At present, the understand- 
ing of power corrections has a solid theoretical 
status only for the special class of processes like 
e + e~ — > hadrons and deeply inelastic scattering 
(DIS) for which the analysis based on the opera- 
tor product expansion (OPE) is applicable. The 
OPE fixes the structure of power corrections (as 
l/Q 4 for the total cross-section of e + e~ annihila- 
tion and l/Q 2 for the structure functions of DIS) 
and allows to identify the corresponding scales A p 
as matrix elements of higher twist composite op- 
erators in QCD. The understanding of power cor- 
rections in the processes which do not admit the 
OPE is. the subject of intensive discussions (see 
review 111 and references therein). 

An important example of the process, which 
does not admit the OPE is the Drell-Yan (DY) 
process Ha + hs — > + X. Here, the lepton 

pair with invariant mass Q 2 is created in the par- 
tonic subprocess q + q — > /-t + /i~(Q 2 ) + X of anni- 
hilation of two quarks with invariant energy s and 
we choose the kinematics such that r = Q 2 / s ~ 1. 
One of the reasons for this is that as r — > 1, the 
final state in the partonic subprocess consists of 
the lepton pair with the energy Q accompanied 
by a soft gluon radiation with the total energy 
(1 — t)Q/2. Due to the presence of two different 
scales the partonic cross-section gets large pertur- 
bative (Sudakov logs in ln(l — r)) and nonpertur- 
bative (as inverse powers-of (1 — t)Q) corrections 
induced by soft gluons Em. Calculating the mo- 
ments <tn ~ Jq dTT N da m I dQ 2 at large N ~ 



and subtracting collinear divergences in the DIS- 
scheme one can represent the leading term of the 
expansion of <7jv in powers of l/Q as 

lna N = ^a s fc (Q/A0C fe (lnA0 + (A VY N/Q) P , 

k 

(1) 

where the coupling constant is defined at the char- 
acteristic soft gluon scale, Ck contains the powers 
of Sudakov logs plus finite, In. JV, contributions 
and the second term represents the leading power 
correction. As N increases, the energy of soft glu- 
ons in the final state decreases as Q/N toward the 
infrared region in which perturbative expressions 
should fail. Indeed, it is well-known that the per- 
turbative expansion in (Q) is not well-defined due 
to factorial growth of coefficients Ck ~ Pgk\p~ k 
at higher orders in a s . This makes the pertur- 
bative series non Borcl summable and induces an 
IR renormalon ambiguity into perturbative con- 
tribution to (|l|) at the level of power correc- 
tions, 5 m a N ~ p oa J Q / N ) ) = ( a QCdN/Q) p . 
Then, IR renormalon ambiguity of perturbative 
expression is compensated in ([j]) by ambiguity in 
the definition of the scale A DY . Thus, examing 
the large order behaviour of the perturbative series 
one can determine the level p at which the leading 
power correction appear in (Q). However, since it 
is impossible to calculate the coefficients Ck ex- 
actly to higher orders in a s one has to find a rea- 
sonable approximation to Ck which would allow to 
identify the k\ factorial growth of the coefficients. 
To this end different schemes were proposed ouB. 
Being applied to the Drell-Yan cross-section they 
predict that there are no N/Q corrections to (0) 
and the leading renormalon contribution has the 
following form 

6 m a N = 0-(N/Q) + (A OY N/Q) 2 . (2) 

We would like to stress that these schemes are ap- 
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proximate and although they were tested in the 
case of c + c~ annihilation and DIS, it is not clear 
yet whether they are applicable in the Drell-Yan 
process and whether they really predict the lead- 
ing power correction to the cross-section. 

2 Wilson line approach 

A different approach to the analysis of power cor- 
rections has been proposed in Ref.cl. It is based on 
the remarkable property of soft gluons that their 
contribution to the Drell-Yan cross-section can be 
factorized with a power occuracy into universal 
factor having the form of a Wilson line. Similar 
to the OPE approach, the power corrections can 
be identified with the matrix elements of certain 
operators built from the Wilson line. 

The Wilson line operator appears in the Drell- 
Yan partonic subprocess as an eikonal phase of 
quarks interacting with soft gluon radiation. Com- 
bining the eikonal phases of quark and antiquark 
together one gets the expression for partonic cross- 
section as 



dQ 2 



= \{n\TU OY (0)\0)\ 2 S(E n (1 - r)Q/2) 



where summation is performed over all possible fi- 
nal states consisting of an arbitrary number, n, of 
soft gluons with the total energy E n = (1 — t)Q/2 
and T stands for time-ordering of gauge fields. 
Here, U DY (0) = Pexp(i J c ^dx ■ A(x)) is the 
eikonal phase of quark + antiquark and A fi (x) 
is the gauge field operator describing soft gluons. 
The contour C(0) in Minkowski space-time cor- 
responds to the classical trajectory of quark and 
antiquark. It goes from — oo to the annihilation 
point along the quark momentum and then re- 
turns to — oo in the direction opposite to the an- 
tiquark momentum. Performing summation over 
final states in the last relation one can obtain the 
representation for the Drell-Yan cross section for 
t — > 1 as a Equrier transformed Wilson line expec- 
tation value El 



da T 



dQ 2 
W(y Q ) 



1 

2 



2tt 



iyo(l-r)Q/2 



TUl(0)TU m (y ), 



(0\W(y Q )\0) 
(3) 



where U DY (yo) is defined at the point y — (yo,0). 
Calculating the moments of this expression one 



obtains the following relation ^ 

o-n = ^tr(0|VK(y o )|0), for y = iN/Q , (4) 

which is valid up to corrections vanishing as TV — > 
oo and which takes into account all Sudakov logs 
and power corrections in N/Q. The scale Q/N en- 
ters into (^) just as a parameter of the integration 
contour in U DY (yo). Eq.([I]) states that the asymp- 
totic expansion of the moments of the Drell-Yan 
cross section is related to the behaviour of the Wil- 
son line expectation value at small yo = 

Perturbative expansion of (^) involves stan- 
dard "short-distance" logarithms a™ln m (|?/o|Ai), 
which are transformed into Sudakov logs for yo = 
iN/Q. The factorization scale fi w Q has a mean- 
ing of the maximal energy of soft gluons contribut- 
ing to the partonic cross-section and after factor- 
ization of soft gluon emission it appears in (|3[) and 
(|) as a ultraviolet renormalization scale of the 
Wilson line operators. As a result, all perturbative 
Sudakov corrections to the Drell-Yan cross-section 
can be effectively resummed using the remarkable 
renormalization properties of Wilson lines. They 
are summarized by the following RG equation 



= 2r cusp (a s )ln(AVQ) 



dln/i 2 



(5) 



with a s = a s (jj,) and r cusp / D Y being some anoma- 
lous dimensions. The evolution of erjv from /x = 
Q/N up to n = Q generates Sudakov logs while 
the initial condition for o~n at /i = N/Q takes 
into account nonperturbative power corrections in 
N/Q. To understand their structure we perform 
the operator expansion of the Wilson line W(yo) 
in powers of yo 



(W(y )) = (W(0))+y Q (dW(0)) 



(G) 



Substituting this expansion into (|J) we find that 
the power corrections to the cross-section are de- 
scribed by matrix elements of the operators de- 
fined as derivatives of the Wilson lineH. Although 
these operators are nonlocal in QCD, one can rep- 
resent them as local operators in an effective non- 
abelian Bloch-Nordsieck theory in which Wilson 
line appears as a result of integration over effec- 
tive quark fields. Solving the RG equation (||) 

"This identity is an analog of the relation between the 
moments of the structure functions of DIS and matrix ele- 
ments of composite operators in the OPE. 
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and using (^) as a boundary condition we obtain 
that the leading power correction to the Drell-Yan 
cross-section appears asl3 

CT P ower = (i8W(0))N/Q + 0(N 2 /Q 2 ) , (7) 

provided that there is no any hidden symmetry 
which would prohibit the appearence of the linear 
~ yo term in the expansion of the Wilson line (|^) 
and would enforce the matrix element (idW(0)) 
to vanish. 

3 IR renormalon analysis 

One possibility to test the presence of the linear 
term in (^|) is to apply the IR renormalon approach 
and identify the source of linear terms inside per- 
turbative series for the Wilson Jjne. At the lead- 
ing order of renormalon calculusu, one can perform 
the explicit one- loop calculation of (W(yo)) with a 
gluon mass A and identify the power corrections as 
nonanalytiral terms in the small (j/o A) expansion. 
One findscl that a s ln(y A) term is absent in agree- 
ment with the Bloch-Nordsieck cancellation, the 
linear term a s (Xyo) also vanishes and the leading 
term is ~ a s (\yo) 2 . This means, that at the lead- 
ing order the contribution of IR renormalons has 
the form (Q). The same result can be obtained in 
the limit of large number of light quark flavours in 
which one performs 1-loop calculation of (W(yo)) 
with the gluon propagator dressed by a chain of 
quark loops □. We notice that in both cases one 
calculates essentially abelian diagrams containing 
the color factors a s CF and a s CV(a s ./V /)" , respec- 
tively. It is well-known that the contribution of 
abelian diagrams to the Wilson line exponentiates 
to higher orders in a s and one can generalize the 
above statement about the absence of linear ~ y 
term to a much larger class of abelian diagrams be- 
yond the leading order of the renormalon calculus. 
However, the natural question rises whether the 
abelian diagrams provide a meaningful approxima- 
tion to the exact expression for the Wilson line or 
may be there is the leading contribution to W(yo) 
coming only from nonabelian diagrams. To give 
an example where the second possibility is realized 
we mention the breakdown of the Bloch-Nordsieck 
cancellation of IR logs in the Drell-Yan cross- 
section at twist 1/Q 4 . The IR divergent part of the 
cross-section has the formH ~ o^CaCf Q _4 lnA 2 
with Ca(Cf) being gluon (quark) Casimir opera- 
tors. It does not appear at order a s and at order 



a 2 it comes only from diagrams with nonabelian 
color structure. One might expect that similar 
phenomenon may happen when one considers lin- 
ear in yo contributions to the Wilson line beyond 
the leading order. 

The straighforward way to check this possibil- 
ity would be to perform 2-loop calculation of the 
Wilson line with a power occuracy in A. However, 
even without going through this complicated cal- 
culation one can employ the operator methods of 
the QCD asymptotic dynamics in order to anal- 
yse the properties of the matrix element 

(dW(0)) = ±-tr{TUl v (0\A)dTU OY (0\A)), (8) 

parameterizing the leading N/Q power correc- 
tion to the Drell-Yan cross-section. The opera- 
tor U DV (yo\A) describes the eikonal phase of fast 
quark and antiquark annihilating at the space- 
time point yo and T— ordering of gauge field op- 
erators A 1 ^ (x) is needed in (||) and (||) to describe 
properly the possibility for quark and antiquark 
to exchange by a virtual soft gluon in the initial 
state. As ajjesult, due to effects of the initial state 
interactional 

TU DY (y \A) ^ U m -(y \A) . 

This does not happen however for single-parton 
initiated process like DIS and it is this phe- 
nomenon that makes different the properties of 
power corrections in DY and DIS. Indeed, if we 
omit the T— ordering in (|J) then it is easy to show 
by direct calculation that the derivative W dU„ v 

" DY D i 

is proportional to the generators of the SU(N C ) 
group and it does not contribute to (||). To 
work out the effects of quark-antiquark correla- 
tions in (W(yo)) one defines two gauge field op- 
erators A^(x) = ^(A fl (x) ± x)) and no- 
tices □ that the field At describes the Coulomb 
gluons while the field A~ is associated with ra- 
diative gluons forming quark and antiquark co- 
herent states. Since [A^' ± (x),Al< ± (y)] = and 
[A^ (x) , A^ (y)] ^ 0, the expansion of the ma- 
trix element (n\TU DV (y \A + + A~)\0) takes the 
following forum 

(n\TU&(y o \A + + A-)\0) = [e i *°] i j jl (9) 

x{n\Ul/(y Q \A-)V:;^(y \A+,A-)\0), 

where the flow of quark color is indicated explic- 
itly. Here, in the r.h.s. the first factor is the non- 
abelian unitary Coulomb phase matrix and the 
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cikonal phase C/ DY depends only on commutative 
operators A~. The residual factor V is unitary, 
V^V = 1, and it describes the correlations be- 
tween Coulomb and radiative gluons. The pertur- 
bative expansion of V looks likeEl 

V = t+ig 3 f abc t a ®t b Jd i xG(x 1 y )A c ll +{x)+O{g i ) 

where effective nonlocal coupling G(x, yo) resem- 
bles the Fadin-Kuraev-Lipatov vertex. In abelian 
theory V — 1 to all orders, but in QCD it gets 
nonabelian corrections starting from order g 3 in 
the strong coupling constant. 

Substituting (0) into (||) one observes that 
nonabelian Coulomb phase is canceled. Moreover, 
if we neglect g 3 — corrections to V, then due to 
unitarity and commutativity of U m (yo\A~), the 
operator in dU DV is again proportional to the 
SU(N C ) generators and its matrix element van- 
ishes identically in Eq.(||). Thus, a nonzero contri- 
bution to (||) comes only from higher order terms 
in V. It immediately follows from the explicit form 
of V that the corresponding Feynman diagrams 
should have a nonabelian color structure and they 
arise starting from order. The set of relevant 
2-loop Feynman diagrams was identified in Ref. u. 
This also explains why linear term did not appear 
in 1-loop calculation of the Wilson line performed 
in Ref.Q. 

The fact that the contribution of the residue 
factor V to the Wilson line W(yo) is nonzeroi-was 
demonstrated by a direct calculation in Ref. u, in 
which a small gluon mass A was used as an IR cut- 
off and only logarithmic terms, ~ ln fc A, were kept. 
It was shown that the leading term W(0) of the ex- 
pansion of the Wilson line in Eq.(^) gets a higher 
twist correction ~ o^CpCAQ^hiX, which leads 
to the breakdown of Bloch-Nordsieck theorem for 
the Drell-Yan cross-section but is nevertheless con- 
sistent with the Kinoshita-Lee-Nauenberg (KLN) 
theorem. To find the linear term in (^|) one has 
to perform the calculation of the same diagrams 
with a power accuracy in A, although interpreta- 
tion of the linear term ~ A to order as a "true" 
power correction is not obvious. One may dress 
instead the gluon propagators by a single chain of 
quark loops, resum all 2-loop diagrams of order 
a^CFCA(ct s Nf) k in the large Nf limit and iden- 
tify the contribution of the leading IR renormalon. 



4 Conclusions 

We conclude that the effects of the initial state 
interaction, that is correlations between Coulomb 
and radiative gluons, provide a new source of N / Q 
power corrections in the Drell-Yan process. In the 
Wilson line approach they are described by the 
matrix element of the derivative of the Wilson line 
operator in Eq.fQ). In the IR renormalon analy- 
sis these effects could manifest themselves only in 
nonabelian Feynman diagrams at higher orders in 
a s . The presence of N/Q corrections in the Drell- 
Yan cross-section is consistent with the KLN the- 
orem. Summation over degenerate initial partonic 
states in (||) and (||) leads to the expression for the 
KLN cross-section, cr* LN = J2k( k \ W ( iN /Q)\ k )> 
where sum goes over an arbitrary number x>f in- 
coming soft gluons. As was shown in Ref. □, the 
initial state interaction effects can be analysed in 
cr^ r LN on the same footing as the final state effects 
thus resulting in the cancellation of N/Q correc- 
tions in cr^ LN . 
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